Surface-Enhanced Plasmon Splitting in a Liquid-Crystal-Coated Gold Nanoparticle 
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We show that, when a gold nanoparticle is coated by a thin layer of nematic liquid crystal, the 
deformation produced by the nanoparticle surface can enhance the splitting of the nanoparticle 
surface plasmon. We consider three plausible liquid crystal director configurations in zero electric 
field: boojum pair (north-south pole configuration), baseball (tetrahedral) , and homogeneous. From 
a calculation using the Discrete Dipole Approximation, we find that the surface plasmon splitting 
is largest for the boojum pair, intermediate for the homogeneous, and smallest for the baseball 
configuration. The boojum pair results are in good agreement with experiment. We conclude that 
the nanoparticle surface has a strong effect on the director orientation, but, surprisingly, that this 
deformation can actually enhance the surface plasmon splitting. 
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Colloidal suspensions and emulsions in anisotropic host 
media, or coated with anisotropic media such as liquid 
crystals have attracted considerable recent interest, since 
these systems exhibit a novel class of colloidal interac- 
tions [1| or the possibility of forming tetravalent binding 
sites of colloids [2| . The alignment of the principal axis of 
anisotropic media (known as the director) around a col- 
loidal surface is a important feature; this new colloidal 
interaction comes from homeotropic alignment, perpen- 
dicular to the surface, while tetravalent binding sites from 
tangential alignment to the surface. 

These systems have also been studied for optical appli- 
cations USIEQ- For example, related to the emerging 
nanotcchnology, a recent experiment has demonstrated 
nanoscale control of optical properties on a length scale 
much less than a light wavelength [jj. The authors ob- 
served the surface plasmon splitting of a gold nanopar- 
ticle coated with a thin layer of nematic liquid crystal: 
the surface plasmon resonance frequency in the scattering 
cross section was shown to depend on the angle between 
the director and the polarization of incident light. They 
also showed that, by rotating the director using a static 
electric field, they could control the change in surface 
plasmon frequency. 

A recent calculation showed that this method could 
indeed produce a measurable change in the surface plas- 
mon frequency, assuming that the director was oriented 
in the same direction everywhere within the coating or 
the host media. But in reality, the director field should be 
influenced near the metal surface , and thus the align- 
ment of the liquid crystal director field near the nanopar- 
ticle surface may reduce the surface plasmon splitting to 
hinder some possible applications of this system which 
involve electrical control of its optical properties. Thus, 
it is important to investigate quantitatively how an inho- 
mogeneous director field may influence the optical prop- 
erties of these coated metal particles. 

In this Letter, we report on novel optical phenomena 
arising from the surface deformation (or alignment) of 
the liquid crystal director field. Our main finding is that 
this deformation can actually enhance this splitting of 



the surface plasmon frequency. We consider three plau- 
sible morphologies for this director field, and, using the 
Discrete Dipole Approximation (DDA) [3, |9(, compute 
the deformation-induced splitting of the surface plasmon 
frequency into two different frequencies polarized paral- 
lel and perpendicular to the director. For the expected 
"boojum pair" configuration, the splitting is enhanced, 
relative to that of a uniform director field, and agrees 
very well with the measured splitting Potential im- 
portance of this work will be discussed below. 

First, we discuss several possible morphologies for the 
liquid crystal director field in an nematic liquid crystal 
layer on a gold nanoparticle. The actual morphology will 
depend strongly on the exact experimental setup. We 
simply assume that, as illustrated in Fig. 1(a), a spher- 
ical coated gold nanoparticle in an air host sits on the 
substrate which is coated with a thin alignment layer, 
so that near the substrate the director will be aligned in 
the x direction. We also assume that the liquid crystal 
director is aligned parallel to the surface at both inter- 
faces. This assumption is reasonable: for a metal-liquid 
crystal boundary, a suitable surface treatment can mod- 
ify the boundary condition to the tangential boundary 
condition (director parallel to the surface) Q, |n| . Also if 
the coating is thin, we can assume a tangential boundary 
condition at the liquid crystal-air interface ^3 • Since the 
anchoring force at the nanoparticle surface can be mod- 
ified by suitable surface preparation, the director may 
sometimes satisfy a more general boundary condition, 
which will not, however, be considered in this Letter. 

The standard analytical description of liquid crystal 
morphology is based on the Frank elastic free energy [lOj . 
If there is no external electric field, assuming the tangen- 
tial boundary condition there are two likely morphologies 
for a spherical nematic liquid crystal layer on a spheri- 
cal nanoparticle. For a relatively thick layer, the director 
field is expected to form a north-south pole type (or boo- 
jum pair) morphology. But if the nematic liquid crystal 
layer is very thin, so that it can be treated like a two- 
dimensional surface, the director field is expected to take 
a baseball- like morphology Q. These two possibilities 
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are shown in Fig. 1. We will calculate the scattering 
cross section for an liquid-crystal-coated gold nanoparti- 
cle, considering both possibilities. 

In our actual optical calculation, we have simplified 
these two morphologies for computational convenience. 
For the boojum pair morphology, we assume that the 
director field is everywhere parallel to the surface and 
points along a latitude between north and south pole. 
This assumption is reasonable because of anchoring, 
which forces the director field to be everywhere parallel 
to the particle surface. The liquid crystal director orien- 
tation is also affected by the alignment layer in the sub- 
strate. If we simply assume a tangential boundary condi- 
tion where the liquid-crystal-coated nanoparticle touches 
the substrate, this condition limits the pole positions. For 
example, in the boojum pair configuration, the alignment 
layer forces the line between north and south poles to be 
oriented in the x direction, i. e. parallel to the substrate. 

For the baseball morphology, there arc four poles, at 
the vertices of a tetrahedron. Because of the substrate- 
imposed anchoring, we assume that two of these are lo- 
cated in lower hemisphere, and that the line between 
these two poles is parallel to the x axis, while the other 
two are located in upper hemisphere, joined by a line par- 
allel to the y axis. To use the DDA, we further simplify 
this configuration by assuming that the director fields in 
the lower and upper hemispheres are both homogeneous, 
and point in the x and y directions, i. e. parallel to the 
substrate surface in both cases. 

To calculate the scattering cross section of the liquid- 
crystal-coated nanoparticle, we use an approach known 
as the DDA [a, |2( • Given the director fields of the liquid 
crystal coating, we calculate the scattering cross section 
Cscat of a single nanoparticle as if none of the other par- 
ticles were present. This amounts to neglecting multiple 
scattering corrections among the different nanoparticles, 
consistent with the experiment of Ref. ||. Also, we ig- 
nore the contribution of the substrate to the scattering 
cross section (although we still include its effects in align- 
ing the liquid crystal director around the nanoparticle as 
described above). Thus, C sca t can be calculated without 
an orientational average over the scattering object. 

As originally formulated, the DDA allows calculation 
of properties such as the extinction or scattering coeffi- 
cient of an irregularly shaped object having a complex, 
frequency-dependent dielectric constant e(w), embedded 
in a homogeneous medium of real dielectric constant e^. 
In this approximation, the object is taken as an array of 
point dipoles (i = 1, . . . , N) with dipole moments P; and 
polarizabilities on, located at positions rj. The extension 
of the DDA to treat an irregularly shaped object which 
has an anisotropic and spatially inhomogcncous dielec- 
tric tensor can be easily done by introducing a position- 
dependent tensor polarizability <Sj, and thus the DDA 
can be applied to our case of an liquid-crystal-coated 
nanoparticle with a spatially varying director field. Thus 



Pi is expressed as 
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the electric field at r^, is due to the incident 
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In the DDA, the product Ajj • P., can be expressed as 
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uj/c = 2ir/\, c and A being 



being the speed of light and wavelength in vacuum. Eqs. 
QJ-© form a coupled set of 3N equations, which can be 
solved for the N dipole moments Pi using the complex- 
conjugate gradient method combined with fast Fourier 
transforms ||. Given the Pi's, C sca t for a given coated 
particle is obtained from thc relation 
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where E se if,i = Ej - E inCii . 

Now we discuss the relation between the local polariz- 
ability tensor a of these point objects and the local di- 
electric tensor e(ui). To carry out the above calculation, 
one needs an expression for the polarizability tensor a. 
For a material with an isotropic scalar dielectric function 
e(<j), one simple means of connecting a scalar polarizabil- 
ity a to e is the Clausius-Mossotti relation. This relation 
can be easily generalized to an anisotropic medium by 
replacing the scalar dielectric function and other scalar 
variables by a dielectric tensor and other tensor vari- 
ables 0| . Thus, we can generalize the lattice disper- 
sion relation |(|, used to connect a to e in the isotropic 
case, to the anisotropic case. To obtain the position- 
dependent anisotropic polarizability tensor a, we insert 
the position-dependent anisotropic dielectric tensor into 
this generalized lattice dispersion relation. 

To apply this approach to the liquid crystal morpholo- 
gies discussed above, we need to transform the local di- 
electric tensor from a frame of reference in which it is 
diagonal (the "diagonal frame") into the lab frame, us- 
ing a suitable similarity transformation. In practice, this 
transformation is not difficult, since the local e and a 
are simultaneously diagonalizable. In practice, we find 
that to calculate a given diagonal element of the lattice 
dispersion relation polarizability, we only need the same 
diagonal element of the dielectric tensor. Thus, the diag- 
onal polarizability tensor using lattice dispersion relation 
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can be calculated in this diagonal frame. Also, the lattice 
dispersion relation is invariant under the above similarity 
transformation. Thus, if the dielectric tensor is initially 
diagonal in some frame of reference, we can obtain the 
diagonal polarizability using lattice dispersion relation 
in this frame first, then carry out a suitable similarity 
transformation on this diagonal polarizability to get the 
polarizability in the lab frame. 

To check the convergence of our DDA calculation, we 
considered an liquid-crystal-coated nanoparticle with a 
gold radius of 30nm, surrounded by an nematic liquid 
crystal of 30 nm thickness. We carried out the DDA 
calculations with a variable number of point objects 
representing the liquid-crystal-coated gold nanoparticle. 
This was accomplished by varying the "linear mesh size" 
M = 2R/d, where d is the distance between adjacent 
point objects on a simple cubic mesh, and R is the radius 
of a metal nanoparticle, i. e. 30nm in this case. In Fig. 
2(a), we plot the scattering cross section as a function of 
this linear mesh size. As can be seen, the peak position 
converges well with increasing linear mesh size M, but its 
height converges less well. We observe the same behavior 
in the case of a nanoparticle with Drude dielectric con- 
stant without an liquid crystal coating layer embedded in 
a homogeneous medium as shown in Fig 2(b). Compared 
with the exact Mie theory, the peak position converges 
well even in case of M = 20, while its height converges 
more slowly (within about 5 % in case of M — 120). In 
Fig. 2, we use Drude dielectric function with ujpT = 7.7. 
Thus, our DDA results for the position of the surface plas- 
mon peak from an liquid-crystal-coated nanoparticle are 
robust against changes of the linear mesh size, even if we 
consider scattering from position-dependent anisotropic 
media. Hence, we can calculate the positions of the scat- 
tering peaks with and without an liquid crystal coating, 
and for various coating thicknesses, using the DDA, and 
compare the results to measurements of Ref . [f| . 

As a first calculation with the actual gold dielectric 
function, we studied how the surface plasmon peak po- 
sition changes as a function of liquid crystal layer thick- 
ness, under the assumption that the director field is uni- 
formly oriented in one direction parallel to the substrate. 
Here for the dielectric function of a gold nanoparticle, 
we used a finite-particle-size corrected dielectric function, 
which we can estimate from the tabulated values of bulk 
gold Q. In this case, we calculated the scattering cross 
section for an incident electric field polarized parallel to 
the liquid crystal director, i. e. = 0° in Fig. 1(b). By 
comparing the position of this calculated surface plasmon 
peak to that measured in Ref. using a similar setup, 
we estimate that the thickness of the liquid crystal layer 
in those measurements is ~ 30 nm, comparable to the 
radius of the gold particle. 

In Fig. 3, we plot the position of the plasmon peak ver- 
sus polarization angle of the incident light, for three liq- 
uid crystal coating morphologies. Two of the morpholo- 



gies are the "boojum pair" and "baseball" configurations 
described above, and the third is a director uniformly ori- 
ented in one direction parallel to the substrate. A mor- 
phology with two singular points of the north-south pole 
type shows the biggest peak shift (about 0.023 eV) with 
change in polarization angle; the magnitude of this shift 
makes the best fit to the value (~ 0.030 eV) observed in 
the experiments. By contrast, the value obtained when 
we assume a uniform director field is about 0.011 eV, 
and that obtained for the baseball-like geometry is the 
smallest among the three different morphologies (~ 0.004 
eV). Thus, the surface deformation of the director field, 
which is caused by the gold nanoparticle surface and by 
the boundaries with host medium or the substrate, can 
strongly influence and even enhance the plasmon split- 
ting of this system. This enhancement suggests that the 
optical properties of such nanoscale systems can be con- 
trolled electrically, by inclusion of liquid crystals. 

The present calculations suggest that the liquid crys- 
tal film in the experiments |]| is probably relatively thick 
(~ 30 nm). If the coating were thinner, the stronger sur- 
face perturbation might lead to a baseball-like configura- 
tion as suggested in Ref. @, and hence, a much smaller 
surface plasmon splitting than seen in the experiments. 
We speculate that there may be a boojum pair-baseball 
crossover as the film gets thinner. If so, monitoring the 
splitting as a function of liquid crystal coating thickness 
would provide an optical signature of such a crossover. 
Moreover, if there is a strong suppression of the plasmon 
splitting, as we find in case of a baseball-like configura- 
tion, this would be an experimental indication that the 
coating thickness is suitable for forming tetravalent func- 
tionalized points around the gold nanoparticlcs, which 
can realize a tetravalent chemistry of colloids Q. 

Finally, we briefly discuss a liquid-crystal-coated gold 
nanoparticle in an applied electric field. If we assume the 
tangential boundary condition, two boojum defects are 
induced on the surface of the nanoparticle Q. In this 
case, the morphology of liquid crystal director near the 
surface of the gold nanoparticle, will be the same as the 
boojum pair case in zero applied electric field, and thus 
we expect that the splitting of the peak in the scattering 
cross section should have the same enhancement by sur- 
face deformation as in zero applied field. The effects of 
an electric field, and more general boundary conditions, 
may need further investigation. 

The present results may have a wide range of ap- 
plications in nanoscience and nanotechnology, because 
they show that optical properties of systems containing 
nanoparticles can be sensitively controlled with the use of 
liquid crystals. The transmission and absorption of such 
materials could be tuned by a dc electric field, which 
will alter the liquid crystal director field, or by control- 
ling the surface interactions between metal and liquid 
crystal. This opens up the possibility of new classes 
of tuneable electro-optic materials based on nanoscale 
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metallic particles combined with liquid crystals. More- 
over, the desirable boojum pair configuration may form 
by self-organization, thereby allowing new types of self- 
assembled nano-optical devices. Finally, a tetrahcdral 
baseball-like director configuration is a potential build- 
ing block for a self-assembled three-dimensional photonic 
crystal with a diamond structure. The measured surface 
plasmon splitting may be a valuable diagnostic to deter- 
mine when this baseball configuration forms. 
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FIG. 1: (a) Schematic of the experimental geometry con- 
sidered in these calculations. We assume that a spherically 
coated nanoparticle sits on a substrate. The substrate has 
a thin coating which aligns the liquid crystal director at the 
point of contact in the x direction. Monochromatic linearly 
polarized light (EM) is assumed normally incident as shown. 
The three possible liquid crystal morphologies which we con- 
sider are shown schematically as (b) homogeneous; (c) "boo- 
jum pair"; and (d) "baseball." We calculate the position of 
the surface plasmon peak as a function of the polarization 
angle 6 indicated in Fig. 1(b). 
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FIG. 2: Calculated scattering cross section versus wavelength 
A for (a) an liquid-crystal-coated metal nanoparticle at angle 
9 — 0° for boojum pair configuration with radius 30nm, liquid 
crystal coating thickness 30nm, and varying the DDA mesh 
size, as indicated in the legend; (b) same for a bare metal 
nanoparticle and a radius 30nm in an air host with various 
mesh sizes. Solid line indicates the exact Mie result. Here we 
used Drude dielectric function, for simplicity. 
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FIG. 3: Frequency of peak in scattering cross section for an 
liquid-crystal-coated gold nanoparticle (with the actual gold 
dielectric function with finite-particle-size correction) versus 
angle 6, for boojum pair, homogeneous, and baseball config- 
urations. 



